function density = PDF_cos(z)
density = 1/2*(1+cos(z*pi));
end

